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Abstract. 

Recently, Serfling and Xiao (2007) extended the L-moment theory (Hosking, 1990) to the multi- 
variate setting. In the present paper, we focus on the two-dimension random vectors to establish 
a link between the bivariate L-moments (BLM) and the underlying bivariate copula functions. 
This connection provides a new estimate of dependence parameters of bivariate statistical data. 
Consistency and asymptotic normality of the proposed estimator are established. Extensive simu- 
lation study is carried out to compare estimators based on the BLM, the maximum likelihood, the 
minimum distance and rank approximate Z-estimation. The obtained results show that, when the 
sample size increases, BLM-based estimation performs better as far as the bias and computation 
time are concerned. Moreover, the root mean squared error (RMSE) is quite reasonable and less 
sensitive in general to outliers than those of the above cited methods. Further, we expect that 
the BLM method is an easy-to-use tool for the estimation of multiparameter copula models. 
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1. Introduction and motivation 

The copula method is a tool to construct multivariate distributions and describe the dependence 
structure in multivariate data sets (e.g., Joe, 1997 or Nelsen, 2006). Modelling dependence struc- 
tures by copulas is a topic of current research and of recent use in several areas, such as financial 
assessments (e.g., Malevergne and Sornette, 2003), insurance (e.g.. Drees and Miiller, 2008) and 
hydrology (e.g., Dupuis, 2007). For the sake of simplicity, throughout the paper, we restrict 
ourself to the two-dimensional case. Let [X^^\ X^'^^) be a bivariate random variable with joint 
distribution function (df) 

F(X1,X2) =P(X« <Xi,X(2) <X2) , (X1,X2) GM^ 

and marg inal df Fj (xj) = P (X(j) < Xj) for Xj £ M and j = 1, 2. If not stated otherwise, we 
assume that the Fj are continuous functions. According to Sklar's theorem (Sklar, 1959) there 
exists a unique copula C : — )■ I, with I = [0, 1] , such that 

F{xi,X2) = C{Fi{xi),F2{x2)), for {xi,X2) eM.\ 

The copula C is the joint df of the uniform random variables (r.v.'s) Uj = Fj [X^^^) , j = 1,2, 
defined for {ui,U2) G I^, by 

C{ui,U2) = F{F^^{u,),F^\u2)), 

where is the generalized inverse function (or the quantile function) of a df G. 

A parametric copula model arises for (^X^^\ X^^^) when C is unknown but assumed to belong to 
a class C := {Ce, 9 e O} , where O is an open subset of W for some integer r > 1. Statistical 
inference on the dependence parameter 9 is one of the main topics in multivariate statistical 
analysis. Several methods of copula parameter estimation have been developed, including the 
pseudo maximum likelihood (PML), inference of margins, minimum distance and others, see for 
instance Genest et al. (2009). All these methods use constrained some optimization techniques, 
and may require substantial computational time. In this paper, we present a new estimation 
method of 9 based on the bivariate L-moments that serves as alternative in front of computation's 
time issue and produces results leads to reasonable estimation performances. The multivariate 
X-moments have been introduced by Serfling and Xiao (2007) as an extension of the univariate 
L- moments introduced by Hosking (1990). The L-comoments have interpretations similar to the 
classical central moment covariance, coskewness, and cokurtosis that also possess the features of 
the Z-moments. This extension is useful to solve some problems in connection with multivariate 
heavy-tailed distributions and small samples. As mentioned, for instance, in Hosking (1990) and 
recently in Delicado and Goria (2008), the main advantage of L-moments vis-a-vis of classical 
estimation methods (e.g. least squares, moments and maximum likelihood) is their relative slight 
sensitivity to outlying data and their performance in statistical inference with small samples. 

In this paper we establish a functional representation of bivariate i-moments (BLM) by the 
underlying copula function and propose a new estimation method for parametric copula models. 
By considering multiparameter Farlie-Gumbel-Morgenstern (FGM) and Archimedean copulas. 
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simulation studies are carried out to compare the performance of this method with those of the 
PML, minimum distance (MD) and rank approximate Z'-estimation. 

The L-moment approach is of interest for multiparameter distributions. In the case of one parame- 
ter distributions, it is equivalent to the classical method of moments (see e.g. Hosking, 1990)(since 
the first L-moment corresponds to the expected value). Next we see, that this property applies 
also in the case of multivariate setting for multiparameter copulas. 

The rest of the paper is organized as follows. In Section [21 we briefly introduce the univariate 
and bivariate L-moment approaches. We present, in Section [3l functional representations of the 
bivariate X-moments by copula functions and give some examples. A new estimator of copula 
parameter and its asymptotic behavior are given in Sections HI In Section \5\ a simulation study 
evaluates the BLM performance is given. 

2. Bivariate X-moments 

First we begin with a brief introduction on the univariate L-moments. Hosking (1990) introduced 
L- moments as an alternative to the classical central moments = E {Y — fi)'^ determined 
by the df Fy of the underlying r.v. Y. An L-moment is defined as a specific linear combination 
of the expectations of the order statistics Yi:^ < ... < lfc:/c- More precisely, the kth L-moment is 
defined by 

By analogy with the classical moments, the first four L-moments Ai, A2, A3 and A4 measure 
location, scale, skewness and kurtosis features respectively. The L-functional representation of A^ 
is terms of the quantile function Fy^ is given by (see Hosking, 1998): 



/ Fy^ {u) Pk-1 (u) du, (2.1) 
Ji 



k 

where Pk{u) := Y.Pk/u'^, with pk/ = {-If'^^ik + i)]/ [(^^^ ! (A; - ^)!] is the shifted Legendre 
polynomials (SLP). In the sequel, we will make use of the three first SLP 

Po (u) = 1, Pi (n) = 2n - 1, P2 (u) = 6^^ -6u + l. 



A straightforward transformation in (j2.ip using Pq = 1 and the orthogonality of Pk-i leads to a 
representation in terms of covariance, that is 

( E[Y] k = l; 

Xk = { (2.2) 
[ Cov{Y,Pk^^{FY{Y))) A: = 2,3,... 

L-moments may be used as summary statistics for data samples, to identify probability distri- 
butions and fit them to data. A brief description of these methods is given in Hosking (1998). 
L- moments are now widely used in water sciences especially in flood frequency analysis. Recent 
studies include Kjeldsen et al. (2002), Kroll and Vogel (2002), Lim and Lye (2003), Chebana and 
Ouarda (2007) and Chebana et al. (2009). In other recent work, Karvanen et al. (2002) used 
L-moments for fitting distributions in independent component analysis in signal processing, and 



4 



Jones and Balakrishnan (2002) pointed out some relationships between integrals occurring in the 
definition of moments and //-moments. Hosking (2006) showed that, for a wide range of distribu- 
tions, the characterization of a distribution by its L-moments is non-redundant. That is, if one 
L-moment is dropped, the remaining L-moments no longer suffice to determine the entire distri- 
bution. Recently, Serfling and Xiao (2007) extended this approach to the multivariate case, this 
has already begun to be developed and applied in statistical hydrology by Chebana and Ouarda 
(2007) and Chebana et al. (2009). 

Next we present basic notations and definitions of the bivariate L-moments. Let X^^^ and X^"^^ 
be two r.v.'s with finite means, margins Fi and F2 and L- moments sequences a[.^^ and A^'^'', 
respectively. By analogy with the covariance representation ()2.2p for L- moments, and the central 
comoments, Serfling and Xiao (2007) defined the kth L-comoment of X^^^ with respect to X^"^^ 
by the covariance of the couple of r.v.'s X^^^ and Pk^i{F2{X^'^^)), for every /c > 2, as 

Afe[12] = COV (F2 (^^'^))) ■ 

Thus, the kth L-comoment of X^"^^ with respect to X^^^ is defined by 

Afc[2i] = Cov [x^^\Pk-i (Fi {x^'^))) ■ 

If we suppose that F belongs to a parametric family of df's, then the A;th L-comoment Xf^[i2] 
depends on the parameters relies the margins and the dependence structure between X^'^^ and 
X^'^\ Since we focus only on the estimation of copula parameters, then it is convenient to use the 
kth L-comoment of Fi (^'•^^) with respect to X^'^^ instead of Afc[i2], that is 

K[i2] ■■= Cov (Fi {x^^^))) , k = 2,3,... 

Similarly the kth L-comoment of F2 (X^^^) with respect to X^-^^ is given by 

K[2i] = Cov {f2 ,Pfc_i (Fi (^('^))) , A: = 2,3, ... 

For the sake of notation simplicity, we set 

4[i2] = A(fc+i)[i2] and (5fc[2i] = A*^^;^)[2i]' k = 1,2, ... (2.3) 

If the copula C is symmetric in the sense that C {u,v) = C {v,u) , then (5fc[i2] = ^^[21], for each 
k = 1,2... We call the coefficient (5fc[i2] by "the kth bivariate copula L-moment" of X^^^ with 
respect to X^'^\ similarly we call 5fc[2i] by the A;th copula L-moment of X^'^^ with respect to X^^\ 

In application, we will often make use of the three first bivariate copula X-moments, that is: 

2Cov (Fi(X«),F2(X(2))) 
-6C0V (Fi (X(i)) , F2 (X{2)) (1 - F2 (X(2)))) 
Cov (Fi (X(i)) , 20F| (X(2)) - 30F| (x(2)) + 12F2 (x(2)) - l) . 



'^l[12] = 
^2[12] = 
^3112] = 
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3. Bivariate copula representation of kth copula i-moment 

Theorem 13 ■ 1 1 b elow gives a representation of the fcth bivariate L-moment in terms of the underlying 
copula function. This result provides a new estimate of bivariate copula parameters. 

Theorem 3.1. The kth bivariate copula L-moment of X^^^ with respect to X*-^-* may be rewritten, 
for each k > 1, as 

^k[i2] = I ^ (C (til, ^2) - U1U2) duidPk iu2) , (3.4) 

or 



Sk\12] = / UiPk {U2) dC iui,U2) . 

Jp 

Observe that (5i[i2] = 5i[2i] = p/S where p is the Spearman rho p, coefficient defined in term of 
copula C by 

p = 12 uiU2dC (ui, U2) — 3, (3-5) 



P 

(see Nelsen, 2006, page 167). 

In view of Theorem 13.11 according to our needs, we can construct a system of equations that will 
serve to the estimation of multiparameter copula models. For this reason, the proposed estimator 
is more likely to be used for the multiparameter copulas. In the case of the one-parameter copulas, 
it is equivalent to the rho-inversion method (see (j4.22p below). Indeed, suppose that we are dealing 
with the estimation of one dimension parameter of a copula model, then it suffices to use one of 
the kth bivariate copula i-moment, says 5i[i2]- In the case of r— dimension parameters we have to 
take the r first bivariate copula Z-moment, so we obtain a system of r equations with r unknown 
parameters. Then, by replacing the coefficients 5fc[i2], ^ = 1, by their empirical counterparts, 
we obtain estimators of the r parameters. Indeed, suppose that r = 3 and C = Cg, 9 = {9i, 62, ^3) , 
then from Theorem 13. H the first three bivariate copula L-moments of X^^^ with respect to X^^^ 
are 

^i\i2] = 2 [ Ce {ui,U2)duidu2 - \ 
Jp 2 

<^2fi2l = 6 / (2n2 - 1) Ce (^1,^2) duxdu2 - \ 
Jp 2 

<^3[i2l = / (60n2 - 60n2 + 12) Cg (ni, U2) duidu2 - 

Jp 2 

Since, in general the copula of {X^^\ X^"^^) is not the same as of {X^'^\ X^^^), the corresponding 
parameters could be estimated accordingly to (5fc[i2] or 5fc[2i]- 

Next we present applications of Theorem 13.11 to parameter estimation of two popular families of 
copula, namely the FGM and Archimedean copulas. 

3.1. FGM families. One of the most popular parametric family of copulas is the FGM family 
defined for \a\ < 1 by 

Ca {ui,U2) = U1U2 + aUiU2UlU2, 0<Ui,U2<l, (3.6) 
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with Uj := 1 — Uj, j = 1,2. The model is useful for the moderate correlation which occurs 
in engineering and medical applications (see, e.g., Blischke and Prabhaker Murthy, 2000 and 
Chalabian and Dunnington, 1998). The Pearson correlation coefficient p corresponds to the model 
(j3.6|) can never exceed 1/3, (see, e.g., Huang and Kotz, 1984). In order to increase the dependence 
between two random variables obeying the type of FGM distribution, Johnson and Kotz (1977) 
introduced the (r — l)-iterated FGM family with r-dimensional parameter a = (ai, ...,ar) ■ 

r 

Ca {ui,U2) = UiU2 + ^ aj (niU2)^-'''^'^^ (uiU2)^-^^^^^/^^ , 
i=i 

where [z] denotes the greatest integer less than or equal to z. For example, the one-iterated FGM 
family (Huang and Kotz, 1984) is a two-parameter copula model: 

Cai,a2 {ui,U2) = U1U2 {1 -|- aiUiU2 + a2UiU2UlU2} ■ (3.7) 

The range of parameters (ai, 02) is given by the region 

7^ := I (ai, 02) ,\ai\ < l,ai + 02 > —I, 012 < - 3 — ai -|- (9 — 6ai — ^afY^^ j- . (3.8) 
The maximal reached correlation for this family is 

Pfgm = 0.42721, for (ai, 02) = (-1 + 7/^13, 2 - 2/V13) , (3.9) 

and the minimal correlation is p^^M ~ "^/^ (oi) 0^2) = (—1, 0) . The two-iterated FGM family 
is given by 

Cai,a2,a3 {ui,U2) = U1U2 |l "I" aiUlM2 + a2UiU2UlU2 + a3UiU2 (miM2)^| , 

and it has been discussed by Lin (1987). 

According to Theorem 13. H we may give explicit formulas of bivariate copula L-moments for the 
FGM, the one-iterated FGM and the two-iterated FGM. Since the number of parameters equals 
k G {1, ...,r} , then we are dealing with first k bivariate copula L- moments that will provide a 
system of k equations and therefore a tool for the estimation of the parameters of the copulas. 
We have : 

• The first bivariate copula L-moment of FGM family Ca is 

'^i[i2] = a/18 

• The two first bivariate copula X- moments of one- iterated FGM copula Cai,a2 'ire: 

<^i[i2l = ai/18 + 02/72 

(3.10) 

^2[12] = "2/120 

The three first bivariate copula L-moments of two- iterated FGM copula Co,i^a2,a3 sue: 

h[i2] = ai/18 + 02/72 + 03/450 

'^2[12] = 02/120 
'53[12] = -03/1050 
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3.2. Archimedean copula families. The Archimedean copula family is one of important class 
of copula models that contains the Gumbel, Clayton, Frank, ... (see, Table 4.1 in Nelsen, 2006, 
page 116). In the bivariate case, an Archimedean copula is defined by 

C{u, v) = ip~^ {ip{u) + ip{v)) , 

where (/? : I — t- R is a twice differentiable function called the generator, satisfying: (p{l) = 0, 
if' (x) < 0, if" (x) > for any x £ 1/ {0, 1} . The notation ip'''^ stands for the inverse function 
of if. For examples, the three generators ipg (t) = (— In ((1 — 6 (I — t)) /t)) , ip,^ (t) = {t~'^ — 1) /a 
and ipjs (t) = (— Int)'' define, respectively, the one parameter Frank, Clayton and Gumbel copula 
families. For example, the Gumbel family is defined by 

C^{u,v) =exp(^- {-\nuf + {-lnvf ^^^^ , f3>l. (3.11) 

For more flexibility in fitting data, it is better to use the multi-parameters copula models than those 
of one parameter. To have a copula with more one parameter, we use, for instance, the distorted 
copula defined by Cr {u,v) = {C (T (n) ,r (v))) , where F : I — )• I is a continuous, concave 
and strictly increasing function with F (0) = and F (1) = 1. Note that if C is an Archimedean 
copula with generator ip, then Cr is also Archimedean copula with generator generator ipoT. For 
more details see Nelsen (2006), page 96 . As example, suppose that F = F^gi with F^^ (*) = 
exp (t~^^ — l) , /32 > and consider a Gumbel copula C/j^ with generator (p/^i (t) = (— Int)'^^ , 
/3i > 1. Then the copula 0/3^^/3^ {u,v) = F'^ {C^^ (F^^ {v))) given by 

Cp,,p, (u, v):={( (u-^^ - l) ' + (v-^' - l) ) +1 ' (3-12) 



which is a two-parameter Archimedean copula with generator (pi3-i^^i32 (t) '■= {t — l)'^^ . (see 
Nelsen 2006, page 96). 

To have the two first bivariate copula iy-moments correspond to C/j^^/jj, we apply Theorem 13.11 to 
get the following system of equations: 

'^i[i2] = 2 /J Jq (C/3i,/32 (u, v) - uv) dudv, 

I I (3.13) 
^2[i2] = 6 /o /o (2u - 1) (C/3,,/32 - dudv. 

In this case we cannot give explicit formulas, in terms of {'^i[i2] > ^2[i2] } i ^^v the parameters {/3i, /32} , 
however for given values of the bivariate copula i-moments, we can obtain the corresponding values 
of {/3i,/32} by solving the previous system by numerical methods. 

Remark 3.1. The previous system provides estimators for copula parameters by replacing the bi- 
variate copula L-moments by their sample counterparts. This is similar to the method of moments 
(see Section^. 

4. Semi-parametric BLM-based estimation 

The aim of the present section is to provide a semi-parametric estimation for bivariate copula 
parameters on the basis of results of Section [3j Suppose that the underlying copula C belongs to 
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a parametric family Cq with 6* = (6*1, • • • , ^j.) G O an open subset of W satisfies the concordance 
ordering condition of copulas (see, Nelsen, 2006, page 135), that is: 

for every 

^(1) _^ 0(2) ^ ^^^^^ <)Cg(2). (4.14) 

The above inequalities mean that for any {u,v) G I^, Cg(i) {u,v) (> or <) Cg(2) (u, w) . It is clear 
that condition ()4.14p implies the well known identifiability condition of copulas: 

for every 9^'^\e^^^ G O : 6^^^ / 9^^^ =^ Cgm + Cgi2). (4.15) 

Identifiability is a natural and even a necessary condition: if the parameter is note identifiable 
then consistent estimator cannot exist (see, e.g., van der Vaart, 1998, page 62). A wide class of 
copula families satisfy condition (j4.14p , and therefore (|4.15p , including the iterated FGM and the 
archimedean models. Indeed, the iterated FGM family is linear with respect to their parameters, 
then by a little algebra we easily verify this condition. For two-parameter Gumbel family, the 
condition ()4.14p is already checked in Nelsen (2006), page 144, example 4.22. Consider now a 
random sample (x^^^ , X^'^^ j from the bivariate r.v. (X(1),X(2)). For each j = 1,2, let 

V / i=l.n 



^tn '■~ ''^Pj-n/ {n + 1) denotes the rescaled empirical df corresponds to the empirical df 



i=l,n 

■= n / ( -I- 1 I HpTiotf 

j:n 

n 

i=l 

We are now in position to present, in three steps, the semi-parametric BLM-based estimation: 
• Step 1: For each k = 1, ...,r, compute 

n 

= --^E^^n (^f ^) (Xf ))) . (4.16) 

i=l 



given in equation (|2.3p . 

Step 2: Using Theorem 13.11 to generate a system of r equations given by equation 

for k = 1, r. 

Step 3: Solve the system 

<^l[12] {(^l,---,0r) = ?1[12] 
'^2[12] (^1, ^'r) = '^2[12] 



K[12] {(^1, ■■■■.Or) = 5r[12]- 



(4.17) 



The obtained solution Q^^^^ ■= ( 6*1, 9^. I is called a BLM estimator for 9 = {9i, 9r) ■ 



The existence and the convergence of a solution of the previous system are established in Theorem 
KT[ (see Section . 

As an application of the BLM based estimation, we choose the one- iterated FGM copula Cai,a2 
given in ()3.7p and propose estimators for the parameters (01,0:2) noted (01,02) . For this family, 
using (|3.10p . system (|4.17p becomes 

O1/I8 + 02/72 = di[i2], 
02/120 = ?2[12], 
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where S^i2], k = 1,2 are given in (|4.16|) . Therefore 

Si = 18(5i[i2] - 30(52[i2], 

32 = 120?2[12]- 

4.1. BLM as a rank approximate Z-estimation. Tsukahara (2005) introduced a new estima- 
tion method for copula models called the rank approximate Z-estimation (RAZ) that generalizes 
the PML one. The BLM method may be interpreted as a RAZ estimation. Let ^{■■,9) be an 
M''-valued function on I^, called "score function", whose components (•; 9) satisfy the condition 

/ '^j{ui,U2;0)dC{ui,U2) = 0, for j = 1, ...,r. 
Any solution 0^^^ in 9 of the following equation 

n 
1=1 

is called a RAZ estimator. There may not be an exact solution to equation (|4.18p in general, 
so in practice, we should choose 0^^^ to be any value of for which the absolute value of the 
left-hand side of equation (j4.18p is close to zero. It is worth mentioning that if the copula Cq is 
absolutely continuous with density cg, then the function ^ = ce/cg, with cg = {dcg/d0k)j=i 
leads the PML based estimation (see for instance Genest et al, 1995). The existence of a sequence 
of consistent roots of Z-estimation in this context is discussed in Theorem 1 in Tsukahara (2005). 
One of the main question of RAZ-estimation is the choice of the score function ^ producing, in 
a certain sense, the best estimator. In Section ([5]), we show that the the copula L-moment score 
(CLS) functions ()4.19p improve the concordance score functions in terms of bias and root mean 
square error (RMSE). 

The univariate L-moments are a special case of the L-statistics with a specific choice of the weight 
coefficients. This fact makes the L-moments interpretable and popular. In a similar way, the 
BLM-based estimator represents an important special case of the RAZ one. In our case, the 
score function ^ (ui, n2; 9) corresponds to the function given in (j4.19p below. Recall that from 
Theorem 13.11 we have 



h\12](.^)= / UiPk{u2)dCg{ui,U2) , k=l,2..., 

Jp 

and define the CLS functions by 

Lk{ui,U2;9) := uiPkiu2) - 6k[i2]i9) , k = l,...,r, (4.19) 

satisfying Jp {ui,U2] 9) dCg {ui,U2) = 0, k = 1, r. Then the RAZ estimator corresponding to 
the CLS function L = (Li, ...,Lr) is a solution in of the system 

n 

1=1 

that is 

n 

E^i- {^^'^) {^L (^f )) " '^'^'^•[121 = 0, = 1, ...,r, (4.21) 
1=1 
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therefore 5fc[i2] {0) = Sk[i2]^ k = l,...,r, which in fact the system of bivariate copula L-moments 
given in system (|4.17p . 

On the other hand, the measures of concordance produce also a Z-estimation for copula models. 
Indeed, the most popular measures of concordance (see, Nelsen, 2006, page 182) are Kendall's tau 
(r). Spearman's rho (p) , Gini's gamma (7) and Spearman's foot-rule phi (99) , given respectively 

by 



r(^)=4/ Cg{ui,U2)dCg{ui,U2) -1, 

Jp 

p{9) = 12 j uiU2dCg (til, 112) - 3, 
Jp 

7 (6*) = 4 Ce (til, 1 - til) dui - / (til - Cg {ui,ui)) dui, 
Ji Ji 

if {6) = 1-3 I til - ti2| dCg (til, ti2) ■ 

Jp 

It follows that the concordance score (CS) functions associated to r, p, 7 and (p respectively are 



^i{ui,U2;e) 
^2 {ui,U2;9) 
^3 (^^i,ii2;^) 
^4 {ui,u2;9) 



ACg{ui,U2)-T{e), 

12ttiti2 - 3-/>(0), 

ACg (til, 1 - til) - til + Cg (til, til) 

1 — 3 |tii — n2| — 93 (6*) . 



It is clear that /j2 ^'j (tii, ti2; 0) dC^i (tii, ti2) = 0, j = 1,...,4, then whenever the dimension of 
parameters r = 4, the function ^ = (^1, ...,^'4) provides Z-estimators for copula models. If the 
dimension of parameters r < 4, then we may choose any r functions from ^i,...,^'4 to have a 
system of r equations that provides estimators of the r parameters. 

Tsukahara (2005) also discussed the RAZ-estimators based on Kendall's tau (r) and Spearman's 
rho {p) , called r-score and p-score RAZ-estimators. Suppose that r = 1 and let and p„ be 
the sample versions of Kendall's tau (r) and Spearman's rho {p) . By using the same idea as the 
method of moments, the r-inversion Or estimator and the p-inversion Op estimator of are defined 
by 

Or = r-^ (r„) and Op = p-^ {pn) . (4.22) 
In the case when r = 2, we may also estimate = [Oi, O2) by solving the system 

(4.23) 




We call, the obtaining solution of the previous system, by (r, p) —inversion estimator of 0. Suppose 
that we are dealing with the estimation of parameters (01,02) of the one-iterated FGM copula 
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Cai,a2 ill (|3.7p . Then, the associated Kendall's tau (r) and Spearman's rho (p) are 

r (ai, 02) = 2q!i/9 + 02/18 + aia2/450 

(4.24) 

p (ai, 02) = ai/3 + a2/12. 
The (r, /o)-inversion estimator of parameters (ai,a2) is the solution of the system 

T(Si,a2) = Tn 

p{ai,a2) = Pn- 

Similarly, if we consider the FGM family Cai,a2,a:j we have to add 7— score and score to have 
a system of four equations, we omit details. 

4.2. Asymptotic behavior of the BLM estimator. By considering BLM's estimator as a 
RAZ-estimator, a straight application of Theorem 1 in Tsukahara (2005) leads to the consistency 
and asymptotic normality of the considered estimator. Let 9o be the true value of 9 and assume 
that the assumptions [^^.1] — [-f^-3] listed below are fulfilled. 

• [-f^.l] 60 £ O C W is the unique zero of the mapping 6* — )• fpi^ (ui, n2; 0) dCog (ni, 112) that 
is defined from O to M''. 

• [-f^.2] L (•; 0) is differentiable with respect to 9 with the Jacobian matrix denoted by 

"(9Lfc {ui,U2\9)' 



L (^1,^2; I 



d9e 



L (ui,U2; 



IS 



L (ui, U2', 9) is continuous both in (ui, U2) and 9, and the Euclidian norm 
dominated by a dCg-integrable function h (ui, U2) ■ 
• The rxr matrix ^0 '■= Jp L (ui, U2; ^o) dCg^ (^1,^2) is nonsingular. 

Theorem 4.1. Assume that the assumptions [H.l] — [H.3] hold. Then with probability tending 
to one as n — )• 00, there exists a solution ^-^la/ equation l \4-21 ] which converges to 9q. 

Moreover 

with 

J2o ■■= var < L (^i, 6; ^0) + V / M^- (ui, ^2) (1 {^j < uj} - uj) dCe,, (^1,^2) i , (4.25) 

where (^1,^2) is a bivariate r.v. with joint distribution function Cq^, 

Mi(ni,n2) ■■= {Pk{u2)}k=i,r and M2 (ni, 7x2) := {uidP^ {U2) / du2} ^.^^^^ ■ 
The proof of Theorem 14.11 is given in the Appendix. 

Note that the assumptions [H.l] — [H.'i] correspond to j42, ^44 and Ah in Tsukahara (2005). The 
score function L meets the remaining assumptions A\ and 743. 
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4.3. A discussion on Theorem 4.1. Notice that assumption [H.l] is verified for any parametric 
copula Cg satisfying the concordance ordering condition of copulas ()4.14p . Indeed, suppose that 
there exists 9i ^ Oq, such that 

/ Lfc (mi,M2; 6*1) ("1' ""a) = 0, for every /c e {!,..., r} . (4.26) 

Jp 

Recall that ^^[12] (6*0) = Jp uiPk (^2) dCe^ {ui,U2) and ^^[12] (6*1) = Jp uiPk (M2) dCe^ (ni,M2) and, 
from assumption (j4.14p . Cqq{> or <)Cq^. It follows, by monotonicity of the integral, that 
^fc[i2] (^'i) (> or <) 5^12] {9o) , this implies that 6k[i2] (Oi) - 4[i2] (Oo) / 0, for every k £ {1, ...,r} . 
Observe that jp dCe^ {ui,U2) = 1 then ^^[12] (^1) = Jp, 5k[i2] (^1) dCe^^ {ui,U2) , consequently, since 
Lk {ui,U2;9i) = uiPk iu2) - (5fc[i2i (6*1) , that 



Lk{ui,U2]9i)dCea{ui,U2) = I UiPk {U2) dCgg {ui,U2) - / 6 u 12] {9 1) dCg^^ {ui,U2) 

P Jp Jp 

= h[i2] (.9o) - Sk[i2] {9i) / 0, for every k G {!,..., r} , 

which is a contradiction with equation (|4.26p , as sought. Let's now discuss the rest of assumptions. In 
[H.2] , the continuity and the differentiability with respect to 9 and {ui,U2) of L (•; 9) and L (•; 9) 
are lie with that of copula Cg, which are natural assumptions in parametric copula models. Some 
examples on this issue are illustrated in Fredricks et al. (2007). The second part of [H.2\ and 
[i7.3] may be checked for a given copula model. For example, if we consider the FGM family 
(see (|3.6p ) we get Li {ui,U2; a) = ui {2u2 — 1) — a/18 and {ui,U2; a) = 0, for k = 2, 3. ..Then 
dLi {ui,U2; a) /da = —1/18 and dL^ {ui,U2; a) /da = 0, for k = 2, ... Let qq denote the true value 
of parameter a. It is clear that each compound of L is continuous with respect to a and {ui,U2) , 
L (ui, U2; a) = 1/18, which is C^— integrable function and, Aq = jp L {ui,U2; ao) dCao — —1/18 

which is nonsingular matrix, then the assumptions [H.2] and [H.3] are well verified. By a little 
algebra we get to the corresponding value of that is defined in (j4.25p , and by Theorem 14. II we 
get 

V^(a^^^^-ao) ^AA(0,ag/270 + l/5) , as n ^ 00. 

For the one-iterated FGM family (see (|3.7p ). by letting ai = a and 02 = it is readily to verify 
that 

Li {ui,U2; a, /3) = ui {2u2 — 1) — a/18 — /3/72, 
L2 {ui,U2; a, l3) = ui [6ul - 6u2 + l) - /3/120, 

which, obviously, are continuous with respect to (a,/3) and {ui,U2) and Cq^/j— integrable function, 
and 

" -1/18 -1/72 
-1/120 

Let (aO)/3o) denote the true value of parameter (a,/3) and by calculating the elements of the 
matrix 

^0:= / ^'{ui,U2',Oio,(3o) dCao,l3o, 

Jp 



L {ui,U2;a,f3) 
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An 



we get 

^ -1/18 1/72 
-1/120 

which is nonsingular because its determinant equals 1/2160 ^ 0, therefore [H.2\ and [H.?i\ are also 
verified. Then, in view of Theorem I4.H we have 

-BLM \ ' " ~ 



a 



n 



"0 



V 



pBLM 








as n — )• oo, 



where ^ := (^q ^) • After a tedious computation we get 



1 



270 540 3780 5 



ao/3o 



8640 2160 







+ 



/3o , «o/3o 



+ 



8640 2160 



an 



«o/3o 

105 252 



+ 



17/3g 1 
21000 15 



it follows that 



342ag 327ao/3o 263/3^ 624 
35 70 ^ 280 ~5~ 



240ag 107ao/3o 443/3^ 



+ 



+ 



140 



+ 240 



240a2 ^10700^^ 443^^2^^ 



7 



7 



140 



960ag 400ao/3o 



7 



+ 



7 



35 



+ 960 



Finally, we note that assumptions [H.l] — [H.'i] may be also verified for one and two parameters 
copula families given in (j3.1ip and (j3.12p . respectively, but that requires tedious calculations which 
would get us out of the context of the paper. 

5. Simulation study 

To evaluate and compare the performance of BLM's estimator with other estimators a simulation 
study is carried out with r = 2 by considering Cai,a2 (the one iterated FGM family) and Cfi^^p^ 
(the two parameters Gumbel family) given in (|3.10p and (j3.12p respectively. The evaluation of the 
performance is based on the bias and the RMSE defined as follows: 

1 ^ / 1 ^ 2\ 

Bias - ^ E - ^) ' I^MSE = E (^^ - ^) ' (5-27) 

i=l \ 4=1 / 

where 9i is an estimator (from the considered method) of 9 from the zth samples for A'" generated 
samples from the underlying copula. In both parts, we selected = 1000. We compare the BLM 
estimator with the PML, (r, /o) —inversion (14.23P and MD estimators (see the Appendix). The 
procedure outlined in Section @ is repeated for different sample sizes n with n = 30, 50, 100, 500 
to assess the improvement in the bias and RMSE of the estimators with increasing sample size. 
Furthermore, the simulation procedure is repeated for a large set of parameters of the true copulas 
Cai,a2 aiid Cp^^p^. For each sample, we solve systems (I3.10p and (I3.13P to obtain, respectively. 
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the BLM-estimators (ai.j,a2,i) and (^/3i,i, /32,i j of (01,02) and (/3i,/32) for i = 1,...,N, and the 
estimators q^, ^fc for k = l,2 are given by = Eili «A:,j and = ^ Xlili 

5.1. Performance of the BLM-based estimation. We first select parameters, as the true 
values of the parameters, of Gumbel and FGM copula models. The choice of the parameters 
have to be meaningful, in the sense that each couple of parameters assigns a value of one of the 
dependence measure, that is weak, moderate and strong dependence. In other words, if we consider 
Spearman's rho p as a dependence measure, then we should select values for copula parameters 
that correspond to specified values of p by using equation (j3.5p . Recall that for the FGM family 
Cai,a2^ the dependence reaches the maximum pf^^j^ = 0.42721 in ai = — 1 + 7/^/13 ^ 0.941 
and a2 = 2 - 2/^/l3 w 1.445 (see (f3^ ). So, we may chose (01,02) = (0.941,1.445) as the 
true parameters of FGM family that correspond to the strong dependence. For the true values 
of (0:1,02) corresponding to the weak and the moderate dependence, we proceed as follows. We 
assign a value to the couple (p,oi) such that |ai| < 1, then we solve by numerical methods the 
equation (j3.5p in the region (|3.8p and get the corresponding value to 02- We summarize the results 
in the following table: 



p 


Ol 


02 


0.001 


0.100 





0.208 


0.400 


0.900 


0.427 


0.941 


1.445 



Table 5.1. The true parameters of FGM copula used for the simulation study . 



By the same procedure, we select the true parameters (/3i,/32) of the Gumbel copula Cp^^p^ and 
get: 



p 




h 


0.001 


1 


0.001 


0.500 


1.400 


0.200 


0.900 


2.500 


1 


ters of Gumbe 


copula 



copula used for the simulation study. 



To evaluate the performance of the ELM estimators, we proceed as follows: 

(1) By using the Algorithm in Nelsen 2006 page 41 and Theorem 4.3.7 in Nelsen 2006 page 
129, respectively, we generate twice A'^ samples of size n from each one the considered 
copulas Ca^^cj and Cp^^p^. 

(2) Obtain the ELM estimators (01,02) of (01,02) and ^/3i,/32^ of (Pi, 1^2) ■ 

(3) By computing, for each estimator, the appropriate Bias and RMSE, we compare (Si,q;2) 
and ^/3i,/32^ , respectively, with the true parameters (01,02) and (/3i,/?2) • 

All computations were performed in the R Software version 2.10.1. The results of the simulation 
study are summarized in Tables [7^ and YTM We observed that BLM's method provides, in terms 
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of bias and RMSE, reasonable results, notably when the sample size increases. However, in the 
case of strong dependence for FGM's family when the sample size is small and less than 30, the 
estimation of the first parameter ai is better that of the second one 02- However, for the sample 
sizes greater than 100 the results become reasonable and more better for sample sizes greater than 
n = 500. For Gumbel family the performance of BLM's method looks good even for small samples. 

5.2. Comparative study: BLM, (t, p) —inversion, MD and PML. As the previous Subsec- 
tion, we consider the bivariate two-parameter FGM and Gumbel copula families with the trues 
parameters those given in Tables 15.11 and 15.21 respectively. The simulation study proceeds as 
follows: 

(1) Generate samples of size n = 30, 50, 100, 500 from the copula Cq. 

(2) Assess the performances of the BLM, (r, p) —inversion, MD and PML estimators. 

(3) Compare the BLM, (r, p) —inversion, MD and PML estimators with the true parameter 
by computing, for each estimator, the appropriate criteria given by (j5.27p . 

It is clear, from Tables [73] to l7.10t that the BLM estimate preforms better than the (r, p) —inversion, 
MD and PML ones as far as the Bias is concerned. On the other hand, in the case of small sam- 
ples the (r, p) —inversion, MD and PLM methods give better RMSE than the BLM one. However, 
when the sample size increases, the RMSE of the BLM estimator becomes reasonable. Moreover, 
for the computation time point of view, we observed that the (r, p) —inversion, MD and PLM 
estimates require hours to be obtained, notably when the sample size becomes large, whereas the 
BLM estimate execution time is in terms of minutes. This is a natural conclusion, because the 
(r, p) —inversion, MD and PLM methods use the optimization problem under constraints, while 
the BLM method uses systems of equations. 

5.3. Comparative robustness study: BLM, (r, p) —inversion, MD and PML. In this 
subsection we study the sensitivity to outliers of BLM's estimator and compare with those of the 
(r, p) —inversion, MD and PML ones. We consider an e— contaminated model for two-parameters 
FGM familly by means of a copula from the same family. In other terms, we are dealing with the 
following mixture copula model: 

Cai,a2 {^) '•— (1 ~ ^) Cai,a2 +£C'ai,a2' (5.28) 

where < e < 1 is the amount of contamination. For the implementation of mixtures models to 
the study outliers one refers, for instance, to Barnett and Lewis (1994), page 43. In this context, 
we proceed our study as follows. First, we select (01,02) = (0.4,0.9) corresponds to Spearman's 
Rho p = 0.208 (see Table [5TT]) and chose (aj,a2) = (0)0) to have the contamination model as the 
product copula that is C'q,*^q,* {u, v) = uv. Then we consider four contamination scenarios according 
to e = 5%, 10%, 20%, 30%. For each value e, we generate 1000 samples of size n = 40 from the 
copula Cai,a2 i^) ■ Finally, we compare the BLM, (t, p) —inversion, MD and PML estimators with 
the true parameter (01,02) by computing, for each estimator, the appropriate Bias and RMSE 
and summarize the results in Table 17. Ill We observed that, for example, in 0% contamination 
the (Bias, RMSE) of Si equals (0.044,0.832) , while for 30% contamination is (-0.165,0.835) . We 
may conclude that the RMSE of BLM's estimation in less sensitive (or robust) to outliers, however 
the Bias is not. The same conclusion is for the (r, p) —inversion method but the BLM one is better. 
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For PLM's estimation both the Bias and the RMSE are sensitive, indeed for 0% contamination the 
(Bias, RMSE) of Si equals (-0.238, 0.440) , while for 30% contamination is (-0.328, 0.589) . Both 
the bias and the RMSE of MD's estimation are note sensitive to outliers, then we may conclude 
that is the better among the four estimation methods. However, the computation time cost in 
MD's method is important which is considered as an handicap from practitioners. 

6. Conclusions 

In this paper, a formula of the bivariate Z-moments in terms of copulas is given. This formula 
leads to introduce a new estimation method for bivariate copula parameters, that we called the 
BLM based estimation. The limiting distribution of the estimators given by the BLM method are 
established. Moreover, we compared by simulations the BLM method with the (r, p)-inversion, 
the minimum distance (MD) and the pseudo maximum likelihood (PML) estimators by focusing 
on the Bias and the RMSE. We conclude that the BLM based estimation performs well the Bias 
and reasonably the RMSE. However, BLM's method may be an alternative robust method as far 
as the RMSE is concerned. As final conclusion, it is worth noting that computation's time of the 
proposed method is quite small compared to MD and PML ones. 
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7.1. Proof of Theorem 13.11 The proof is straightforward and based on elementary calculation. 



7. Appendix 



Indeed, since 




is (0, 1) —uniform r.v., then copula's representation of the joint df of the 





then, the covariance of (x^^\ (F2(X(^))) j equals 
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12 



D{u, v) — V 

C (^Fi{u),v^/^^ - v^/^Fi{u)^ dvdu 
{C{u,v)-uv) dv^dF{^{u), 



it follows that 



i=l ^ ^ 



Since 



^fc[12] 



Cov (x^'IpJfJx^'^ 



£=1 ^ ^ 



y^,P£,k {C {u,v) - uv)dv^du 



{C (u, v) — uv) dPk (v) du, 



as sought. 



□ 



7.2. Proof of Theorem 14.11 The existence of a sequence of consistent roots 0^^^^ to (j4.17p or 
(j4.2ip . may be checked by using a similar argument as the proof of Theorem 1 in Tsukahara 
(2005). Indeed, we have only to check the conditions in Theorem A. 10.2 in Bickel et al. (1993). 
Since we are dealing with an asymptotic result, we may consider that, for all large n, without 
loss of generality, that the empirical df Fyn and their rescaled version Ft„ have a same effect. 
Therefore throughout the proof, we will make use of Fj-^ instead of F^^. For convenience we set 

1 " I' 

^n{e) := -J^L (^i:n (^f^) , ^2;™ (^f^) ^o) and $ (0) := J^Liui,U2;e)dCeo {ui,U2) . 

i=l 

By assumption [H.2] , it is clear that the foUowing derivatives exist: 



1>n (0) 



$(0) 



d^n (0) _ 1 

de 

(9) 

de 



1 = 1 
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L{ui,U2;O)dC0Q {ui,U2). 



Next, we verify that 



sup 



(9) - 1> {9) 



'o| < ^ 0, as n ^ oo, 



(7.29) 
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for any real sequence e„ — )• 0. By using the triangular inequality we get 
1 



(9) - (0o) 



< 



n 



E 

i=l 



L(F,,AX^'A,F,,JxP 



L(Fi,„(X«),F2:„(xfM;^o 



Since L is continuous in 0, then 



sup 
therefore 



sup 



ol < r ~^ 0' as n — )• oo. 



(0) - (^o) 

On the other hand, from the law of the large number, we infer that 



'ol <en> = Op (1) 



(7.30) 



as n — )■ oo. 



i=l 



Moreover, in view of the continuity of L in w and since sup^(j) \Fi:n 
almost surely, n — )• oo (Glivenko-Cantelli theorem), we have 



0. 



It follows that 



(^o) - $ (^o) 



0, which together with fOO]) . implies (f7:29|) . 



Conditions (MGO) and (MG3) in Theorem A. 10.2 in Bickel et al. (1993) are trivially satisfied by 
our assumptions [HI] — [H'i] . In view of the general theorem for Z-estimators (see, van der Vaart 

and Wellner, 1996, Th. 3.3.1), it remains to prove that ^/n ~ (^o) converges in law to the 

appropriate limit. But this follows from Proposition 3 in Tsukahara (2005), which achieves the 
proof of Theorem 14.11 □ 

7.3. Minimum distance based estimation. We briefly present the minimum distance (MD) 
base estimation for copula models that possesses a qualitative robustness (Genest and Remillard, 
2008), this will be compared with the BLM method (see Subsection 5.2). Let C be the true 

copula associated to the df of (X^^^X^^^) and suppose that we have a given parametric family 
of copula C := {Cg, 9 £ O} to fit data. Let us define the minimum distance functional T on the 
space of the copula by 

T(C) := argmin/i(C,Ce) . 

Here /x is a distance between probabilities on I^. In the present paper, we consider the Cramer- von 
Mises distance defined by 

/i™(C,C,) := / {C{ui,U2)-Ce{uuU2)fdC{uuU2). 
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Consider now a random sample (x^^^ , X^"^^ ) , from the bivariate random variables (X^^'' , ^^^^ ] 

V / i=l,n 

The joint empirical distribution functions is given by 



1 " 

Fn{xuX2) = <Xi,Xf) <X2] 



n . 



Following Deheuvels (1979), we define the empirical copula by 

Cn (ni,n2) := F„ {¥'1 [ui) ,F~l (us)) , < -ui.-ua < 1. 
The corresponding Cramer-von Mises statistics is 

pP^^(Cn.Ce)= [ {Cn{ui,U2)-Ce{ui,U2)fdCn{ni,U2). 

This may be rewritten into 
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/.^^^ (C„, Co) = n-' U\'^) - Co {p^,uf^ 

i=l 

where U^^^ := (^^i"'^) ; ^ = for each j = 1,2 (see, Genest and Remillard, 2008, eq. 

31). The MD estimator of the parameter 6 is defined by 

e = T{Cn) :=argmin//^^*^ (C„,Ce). 

Note that we may also use the Kolmogorov-Smirnov distance but this is awkward in practice due 
to the supremum norm uses. Also since the Hellinger distance is defined by copula densities, other 
nonparametric estimators of the underling copula are needed (see, Biau and Begkamp, 2005) and 
therefore non-standard computational procedures are required. 

Suppose now that we are dealing with the estimation of parameters of one iterated FGM copula 
family Cai,a2 in (13. 7p by means of the MD method. The MD estimator for a = (ai,a2) noted 
Q^MD pgg^ii^s ]-,y minimizing the function (01,02) — p{Cn,Cai,a2) over the region TZ given in 
()3.8p . Then to solve the previous optimization problem, we will introduce the Lagrange multiplier 
principle, that is we have to rewrite the region TZ into 

n ={(01,02) , £j (01,02) > 0, j = 1,2,3} , 

where ii (oi, 02) := 1 — of, £2 (oi, 02) := oi + 02 + 1 and 



and then minimize the function 



^3(01,02):=- 3 — oi + (9 — 601 — 3oi)^''^ 



02, 



3 

Kn (o,z^) := p(C„,Cq,,,q,2) - '^i^jij (01,02) , 

i=i 

over the whole M^, with o = (oi, 02) G and = (i^i, 1^2, 1^3) £ M'^. So, the new formulation of 
the MD estimator of parameter o is 

(a,v) - 
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We note here that it is difficult, in general, to have an explicit form for a , then only the 
numerical computation can solve this issue. This is observed for the one-iterated FGM family, 
that the optimization problem requires tedious tools. 







P = 


0.001 






P = 


0.208 






P = 


0.427 






ni = 


0.1 




= 


n 1 


= 0.4 


no = 


0.9 


fVi = 


0.911 


no = 


1.445 


n 


Bias 


RMSE 


Bias 


RMSE 


Bias 


RMSE 


Bias 


RMSE 


Bias 


RMSE 


Bias 


RMSE 


30 


0.197 


0.882 


-0.089 


2.417 


0.113 


0.931 


-0.210 


2.804 


0.098 


0.823 


-0.312 


2.861 


50 


0.133 


0.672 


-0.050 


1.781 


0.042 


0.703 


-0.065 


2.074 


-0.051 


0.712 


0.276 


2.197 


100 


0.065 


0.456 


-0.040 


1.105 


0.026 


0.498 


0.048 


1.408 


0.041 


0.513 


-0.055 


1.572 


500 


-0.017 


0.206 


0.041 


0.639 


0.021 


0.215 


0.031 


0.659 


-0.020 


0.308 


-0.031 


0.692 



Table 7.3. Bias and RMSE of BLM's estimator of two-parameters FGM copula. 



p = 0.001 


p = 0.5 


p = 0.9 


/3i = 1 p2 = 0.001 


Pi = 1.4 132 = 0.2 


/3i = 2.5 132 = 1 


n Bias RMSE Bias RMSE 


Bias RMSE Bias RMSE 


Bias RMSE Bias RMSE 


30 0.162 0.994 0.428 1.945 
50 0.134 0.725 0.294 1.107 
100 -0.094 0.697 0.219 0.804 
500 -0.071 0.597 0.107 0.358 


0.214 1.002 0.549 1.421 
0.187 0.695 0.498 0.999 
-0.136 0.619 0.287 0.665 
-0.081 0.489 0.148 0.477 


0.404 0.920 -0.653 1.109 
0.350 0.854 -0.550 0.835 
0.183 0.597 -0.536 0.526 
-0.096 0.395 -0.340 0.480 



Table 7.4. Bias and RMSE of BLM's estimator of two-parameters FGM copula. 



(>l - 0.1 (>2 - 





Bias 


RMSE 


Bias 


RMSE 


Time (h) 


n = 30 


BLM 


0.227 


0.952 


-0.194 


2.882 


0.640 


(r, p)-inversion 


0.458 


1.005 


0.782 


2.157 


0.978 


MD 


0.575 


0.571 


0.494 


0.851 


1.566 


PML 


0.550 


0.552 


0.424 


0.872 


1.033 


n = 50 


BLM 


-0.140 


0.702 


-0.112 


2.193 


1.215 


(r, /9)-inversion 


0.358 


0.958 


0.558 


1.428 


2.856 


MD 


0.468 


0.559 


0.238 


0.846 


3.455 


PML 


0.444 


0.546 


-0.237 


0.840 


2.421 


n = 100 


BLM 


-0.039 


0.565 


0.082 


1.364 


1.847 


(r, p)-inversion 


0.229 


0.664 


0.195 


0.985 


2.548 


MD 


0.125 


0.521 


0.145 


0.684 


6.888 


PML 


0.121 


0.520 


0.131 


0.673 


4.107 


n = 500 


BLM 


0.021 


0.417 


0.071 


0.634 


5.963 


(r, /9)-inversion 


0.084 


0.588 


0.118 


0.748 


11.548 


MD 


0.077 


0.504 


0.086 


0.640 


19.598 


PML 


0.076 


0.502 


0.081 


0.639 


17.073 



Table 7.5. Bias and RMSE of the BLM, (r,p)-inversion, MD and PML estimators 
for two-parameters of FGM copula for weak dependence {p = 0.001). 





oci = 


0.4 


012 = 


0.9 






Bias 


RMSE 


Bias 


RMSE 


Time (h) 


n = 30 


BLM 


0.127 


0.855 


-0.297 


2.668 


0.011 


(r, p)-inversion 


0.102 


1.322 


-0.290 


1.383 


0.516 


MD 


-0.174 


0.777 


-0.322 


1.058 


3.583 


PML 


-0.191 


0.906 


-0.372 


1.261 


0.954 


n = 50 


BLM 


-0.059 


0.755 


0.123 


2.001 


1.035 


(r, /9)-inversion 


0.091 


0.892 


-0.141 


1.272 


2.101 


MD 


-0.173 


0.730 


-0.223 


1.010 


6.428 


PML 


0.122 


0.775 


-0.200 


0.853 


2.823 


n = 100 


BLM 


0.031 


0.715 


0.060 


1.404 


1.920 


(r, /9)-inversion 


0.082 


0.791 


-0.130 


0.942 


2.037 


MD 


-0.130 


0.652 


-0.121 


0.919 


11.217 


PML 


0.090 


0.599 


0.100 


0.794 


3.652 


n = 500 


BLM 


-0.025 


0.300 


0.049 


0.629 


5.205 


(r, /9)-inversion 


0.054 


0.393 


-0.087 


0.701 


8.285 


MD 


-0.071 


0.602 


-0.061 


0.742 


19.210 


PML 


0.047 


0.573 


0.056 


0.632 


16.458 



Table 7.6. Bias and RMSE of the BLM, (r,/9)-inversion, MD and PML estimators 
for two-parameters of FGM copula for moderate dependence (p = 0.208). 









(>2 - 1. 


145 






Bias 


RMSE 


Bias 


RMSE 


Time (h) 


n = 30 


BLM 


0.091 


0.832 


0.402 


2.715 


0.017 


(r, p)-inversion 


0.171 


1.142 


-0.471 


1.229 


0.757 


MD 


-0.142 


0.871 


-0.420 


1.025 


2.083 


PML 


-0.121 


0.927 


-0.415 


1.061 


1.781 


n = 50 


BLM 


0.054 


0.641 


0.300 


1.982 


1.020 


(r, y9)-inversioii 


0.157 


0.997 


-0.321 


1.120 


2.021 


MD 


-0.135 


0.753 


0.351 


0.940 


6.633 


PML 


0.092 


0.892 


-0.307 


1.150 


2.754 


n = 100 


BLM 


0.030 


0.449 


0.090 


1.391 


1.620 


(r, p)-inversion 


0.081 


0.463 


-0.153 


0.931 


3.037 


MD 


0.070 


0.743 


-0.114 


0.904 


9.217 


PML 


0.050 


0.712 


0.102 


0.800 


3.652 


n = 500 


BLM 


0.021 


0.315 


0.046 


0.602 


5.205 


(r, /9)-inversion 


0.071 


0.357 


-0.098 


0.765 


8.285 


MD 


-0.064 


0.541 


-0.054 


0.782 


19.210 


PML 


0.052 


0.472 


0.076 


0.699 


17.458 



Table 7.7. Bias and RMSE of the BLM, (r,p)-inversion, MD and PML estimators 
for two-parameters of FGM copula for strong dependence {p = 0.427). 





/3i 


= 1 


/32 = 


0.001 






Bias 


RMSE 


Bias 


RMSE 


Time (h) 


n = 30 


BLM 


0.174 


0.941 


0.453 


1.854 


1.121 


(r, /9)-inversion 


0.181 


0.782 


0.532 


1.186 


2.021 


MD 


-0.274 


0.546 


-0.698 


1.243 


4.691 


PML 


0.310 


0.335 


-0.593 


0.910 


2.065 


n = 50 


BLM 


-0.157 


0.897 


0.289 


0.977 


1.026 


(t, p)-inversion 


0.184 


0.539 


0.476 


0.629 


2.265 


MD 


0.262 


0.448 


-0.310 


0.759 


5.633 


PML 


0.250 


0.303 


-0.302 


0.815 


2.754 


n = 100 


BLM 


-0.126 


0.530 


0.193 


0.824 


1.920 


(r, /9)-inversion 


-0.177 


0.523 


0.250 


0.619 


2.248 


MD 


-0.161 


0.420 


-0.201 


0.521 


6.285 


PML 


0.151 


0.272 


-0.197 


0.810 


4.153 


n = 500 


BLM 


-0.098 


0.411 


0.114 


0.324 


5.010 


(r, p)-inversion 


-0.235 


0.502 


0.136 


0.503 


7.149 


MD 


-0.181 


0.409 


0.116 


0.376 


14.984 


PML 


0.170 


0.205 


-0.115 


0.619 


13.147 



Table 7.8. Bias and RMSE of the BLM, (r,/9)-inversion, MD and PML estimators 
for two-parameters of Gumbel copula for weak dependence {p = 0.001). 



Ji - J2 - 0.2 





Bias 


RMSE 


Bias 


RMSE 


Time (h) 


n = 30 


BLM 


-0.182 


0.989 


0.593 


1.317 


1.024 


(t, p)-inversion 


0.191 


0.885 


0.655 


1.215 


2.042 


MD 


0.214 


0.985 


-0.525 


1.056 


4.485 


PML 


0.195 


0.524 


-0.423 


1.051 


2.125 


n = 50 


BLM 


-0.134 


0.594 


0.526 


0.994 


1.058 


(r, /9)-inversion 


0.187 


0.512 


0.555 


0.972 


2.212 


MD 


0.181 


0.423 


-0.461 


0.853 


5.588 


PML 


0.177 


0.318 


-0.413 


0.916 


2.859 


n = 100 


BLM 


-0.122 


0.482 


0.272 


0.492 


2.247 


(r, p)-inversioii 


-0.150 


0.421 


0.291 


0.712 


3.153 


MD 


-0.170 


0.439 


-0.269 


0.474 


6.256 


PML 


0.152 


0.293 


-0.275 


0.471 


5.254 


n = 500 


BLM 


-0.101 


0.223 


0.135 


0.312 


6.587 


(r, /9)-inversion 


-0.149 


0.400 


0.221 


0.655 


9.145 


MD 


-0.106 


0.306 


-0.212 


0.355 


14.445 


PML 


0.102 


0.221 


-0.200 


0.317 


13.157 



Table 7.9. Bias and RMSE of the BLM, (r,p)-inversion, MD and PML estimators 
for two-parameters of Gumbel copula for moderate dependence {p = 0.5). 



/3i = 2.5 /32 = 1 





Bias 


RMSE 


Bias 


RMSE 


Time (h) 


n = 30 


BLM 


0.422 


0.954 


-0.740 


1.119 


0.755 


(r, /9)-inversion 


0.786 


1.125 


0.782 


1.175 


1.254 


MD 


0.546 


0.546 


0.592 


0.563 


1.545 


PML 


0.553 


0.551 


0.723 


0.522 


1.765 


n = 50 


BLM 


0.329 


0.817 


-0.635 


0.852 


1.021 


(r, p)-invGrsion 


0.586 


0.983 


0.745 


0.972 


2.045 


MD 


0.321 


0.522 


0.582 


0.552 


2.265 


PML 


0.292 


0.512 


0.551 


0.514 


2.255 


n = 100 


BLM 


0.107 


0.584 


-0.592 


0.713 


1.920 


(r, /9)-inversion 


0.425 


0.812 


0.611 


0.902 


2.153 


MD 


-0.181 


0.501 


-0.578 


0.488 


5.544 


PML 


0.172 


0.482 


-0.545 


0.472 


5.458 


n = 500 


BLM 


-0.066 


0.456 


-0.367 


0.478 


5.205 


(r, y9)-inversioii 


0.123 


0.757 


0.501 


0.694 


9.789 


MD 


0.094 


0.469 


0.408 


0.495 


14.565 


PML 


0.084 


0.465 


0.375 


0.482 


13.425 



Table 7.10. Bias and RMSE of the BLM, (r,p)-inversion, MD and PML estima- 
tors for two parameters of Gumbel copula for strong dependence (p = 0.9). 



(>i - 0, 1 02 - (1.9 



0% contamination 




Bias 


RMSE 


Bias 


RMSE 


BLM 


0.044 


0.832 


-0.141 


2.650 


(r, p)-inversion 


0.053 


0.432 


-0.440 


0.711 


MD 


0.267 


0.270 


-0.456 


0.461 


PML 


-0.238 


0.440 


-0.472 


0.627 


5% contamination 


BLM 


0.046 


0.833 


-0.137 


2.662 


(r, /9)-inversion 


-0.069 


0.431 


-0.479 


0.738 


MD 


0.254 


0.257 


-0.472 


0.475 


PML 


-0.274 


0.407 


0.432 


0.613 


10% contamination 


BLM 


-0.082 


0.811 


-0.155 


2.641 


(t, p)-inversion 


-0.090 


0.393 


-0.461 


0.695 


MD 


0.279 


0.281 


-0.464 


0.468 


PML 


-0.267 


0.506 


-0.429 


0.637 


20% contamination 


BLM 


-0.100 


0.802 


-0.188 


2.585 


(r, /?)-inversion 


-0.130 


0.423 


-0.537 


0.786 


MD 


0.280 


0.282 


-0.472 


0.477 


PML 


-0.268 


0.524 


-0,500 


0,639 




'My/l com. 


;iiiii:iaLioii 






BLM 


-0.165 


0.835 


-0.280 


2.627 


(r, p)-inversion 


-0.179 


0.480 


-0.619 


0.909 


MD 


0.293 


0.266 


-0.458 


0.465 


PML 


-0.328 


0.589 


-0.515 


0.641 



Table 7.11. Bias and RMSE of the BLM, (r,/9)-mversion, MD and PML estima- 
tors for e-contaminated two-parameters of FGM copula by product copula. 



